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The aim of this paper is to show a way of applying Bayesian inference on interval fuzzy data assuming that 
all the time we work with interval probabilities. The proposal is illustrated with an example in Health Care, 
specifically on inferring population annoyance level caused by noise exposure. 


0.1 Introduction 


In this paper we introduce a way of applying Bayesian inference in the presence of linguistic imprecision at 
the same time in evidence, priors and marginals. Linguistic imprecision is represented using interval fuzzy 
and interval probability assignments. Our ideas are illustrated with an example in Health Care, specifically 
on inferring population annoyance level caused by noise exposure. Expressions such as “In the light of these 
results we can say that for a sound level of approximately six or seven, then, with a probability, approximately 
high or perhaps very high, the best alternative is that people living at that place is moderate annoyed because 
of acoustic contamination” are allowed. 


0.2 Imprecise probabilities model 


In the presence of imprecision, instead of a single probability distribution, we can use a class M of priors, 
so that each prior in M is updated by Bayes’ rule, producing a class of posteriors. In the theory of imprecise 
probabilities the lower and upper probabilities of an event or hypothesis H, denoted by P(H) and P(H), are 
defined by P(H) = inf{P(H) | P € M} and P(H) = sup{P(H) | P € M}. This theory denies that any 
precise probability model is reasonable. Statistical conclusions will be expressed in terms of posterior lower 
and upper probabilities or expectations, and these imprecise probabilities should not be regarded as lower and 
upper bounds for some unknown precise probability [i]. 
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Let © be the parameter space and A C O. The lower and upper posteriors of A are obtained minimizing and 
maximizing, respectively, the posterior probability 


a(Alx) =f m*(ebsaa / f «*(opx)a0 (2) 


in the lower/upper class of unnormalized posteriors 7*( |x), which we assume being represented as an interval 
of measures [2]. In order to minimize 7(A|x), we use (Cf. [f]) 


m*(8|x) = (6|x)xoea + U" (9X) Xoce\4 (2) 


in (2), where yp denotes the characteristic function of predicate P, i.e. yp = Lif P is true, and yp = Oif P is 
false. Thus, we assign the minimum possible mass to A and the maximum possible mass to the complement of 
A in O. The lower posterior is 


B(Alx) =f e@xae / [[.reooan [erro (3) 


In order to maximize 7(A|x), we only have to swap /* with u* in (2). 


0.3. First Ideas on Fuzzy Bayesian Inference 


A linguistic variable whose values are words or sentences may be defined by the quadruplet 
(Tname, L(E),E, M), where Lpame is the name of the linguistic variable; L(E) is the term-set or reference- 
set Of Lname, i. e., the finite set of linguistic values that ®name can take, which elements we denote by e;; E€ is the 
universe of discourse or physical domain associated with Ypame; M is a semantic function that associates a fuzzy 
meaning to each linguistic value ez € L(E), i. e., it is an injective mapping from L(E) to F(E), the set of fuzzy 
subsets of €. As it is more and more usual, we denote M (ez) by ez. Remember that the support and the core of 
A € §(€) are supp(A) = {x €E: A(x) > 0} and core(A) = {x EE: A(x) = 1}. 

First works on fuzzy Bayesian inference come from safety project studies in structural reliability researches 
[4]. Given a linguistic value e7, of the evidence, and a set of exhaustive and mutually exclusive hypotheses H; 
(j = 1,...,m), we can compute the likelihood p(ez,|H;) by 


p(ex|H;) = / _, Hex(€)F(€\H)d (4) 


where f(e|H;) is the likelihood density function evaluated at e given the hypothesis H;. We can compute the 
posterior probability by 


p(dlen) = aH plexi) / SO, (w(t )PCen| Ft) (5) 


One problem is the determination of f(e|H,;) (usually approximated as Gaussian or Weibull). On the other 
hand, once f(e|H;) has been determined, another problem is the computation of the integral itself. In order to 
save us from these problems, Yang [5] proposes to deduce f(e|H;) from p(ez|H;): 


f(e|H;) mie p(er|H;); withe =W ca)/ | Le, (e)de 
er EL(E e€supper 


Lastly, we can compute the posterior probability by the Bayes’ rule, 


(dle) = #elHt oH) / SO felt (Hh) (6 


The use of c means that the size W (ez) of the range covered by any linguistic value e,, divided by the area 
of its membership function ji, , is the same for each linguistic value e7,, and is equal to c. 
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0.4 Interval Fuzzy Bayesian Inference (IFBI) 


Our proposal does not depend on any constant and it works with interval fuzzy and interval probability 
assignments. Let e7, be a fuzzy value for the evidence. Let O = {H1,..., H,,} be a complete set of hypotheses, 
ie. a finite set of exhaustive and mutually exclusive hypotheses. We define the likelihood function for any 
continuous value e € € of the evidence (given an hypothesis H;), from the likelihoods of the fuzzy values e;, of 
the evidence (given the same hypothesis H;): 


LUHjle) «> carey ee) P(ez| A) (7) 
We extend (7) to fuzzy values E € §(E) of the evidence by 
L(Hi|E) «Yo ng OE(E)P(er| Hj) (8) 


where ez, : §(E) — §([0, 1]) extends ey : € — [0,1] according to Zadeh’s extension principle (Cf. Footnote [a). 

Let n,m be positive integer numbers. Let I”R be the set of all intervals of dimension n of real numbers. 
Given D C R™, we denote by I” D the set of all intervals of dimension m included in D. Given f : D > R”, 
an inclusion function for f, is any function gf : ID — IR, such that, VJ € I"D, f(J) C of (J). If f 
is monotonic increasing, such as e”, In z, or \/z, then an inclusion function is gf ({ao,a1]) = [f (ao), f(a1)). 
If f is monotonic decreasing, then an inclusion function is gf ({ao,a1]) = [f(a1), f(ao)]. The case of non 
monotonic functions is not difficult if we know the intervals where the function is monotonic. For example, in 
the illustrative example below, we use triangular fuzzy numbers (tfn), because of their ease of computation. 
The membership function of a tfn is characterised by the lower, modal and upper values, i.e. the vertices of the 


triangle: 
p(x; a,b,c) = max (nin (= 1 ——*)..0) (9) 


LEL(E) 


b-—a’’c—b 


The natural extension of a triangular fuzzy number T(z; a, b, c) is 


Teli. be) = Peon, 0.6), T (era, 0, C)| Rares (10) 
+ ata; a, b, C), Eros a, b, c)] Xb<29 
+ [min{T (x9; a, b, €) Tay; a, b, aie lL) Xegcoeer 


Thus, the natural extension of (8) to interval fuzzy values of the evidence is 


L(H,| [E]) « 5> er ([E])p(er| Hy) 


eLrEL(E) 


Le. an interval ( L(H|[El)) L(H;|[E1)) J. If [E] = [Eo, Ei] then ( L(H|[El)) Se ehy) 
and ( L(H;| (E))), - L(H;|F1), calculated according to tt 


The lower and upper likelihoods provides a lower/upper class of inclusion functions for the unnormalized 
posteriors, 


S* = {op" : Vj, al* (Ay [E]) < ob" (Hy| [E]) < ot" (Hjl [E])} (13) 


where /* (H;|[E]) and pu*(H;|[E]) are proportional to ( L(H;| (B])) P(E) and to 
L(H;| (E))) , p(H;), respectively. 


Let S be the class of normalized inclusion functions 4p(H;| |E]) for the posteriors. If we whish a lower/upper 
expression for the class S, we have to minimize and maximize pp(H;| [E]) over S*. To minimize it, we take 


PD (H| (E]) = ol* (Al (El) Xuan, + ot" (Al (E)) xu zn, (12) 


in order to assign minimum possible mass to H; and maximal mass to every H # H;. Thus, we obtain the 
normalized lower posterior, 


i (#,;| [B)) = of* (EE) ie) / ( PGE) +0 
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To maximize p(H;| [E]) simply swap 7* with gi* in (i2}. Swapping p/* with gu* in we obtain the 
expresion for the normalized upper posterior Gu (H;| [E)). 

The following proposition shows that although neither I nor pw are additive over © they both distribute a 
whole mass of two between any subset A C O and its complement O\ A. 


Proposition 1 The lower/upper posteriors satisfy: 


(i) 1 € coregl(O| [E]) N corequ(9| [F]); 
(ii) al (A| [E]) + ol((O\A)| [E]) = one = ot (Al [E]) + o%t((O\ A)| [E]); 
(iii) gl (A| [E]) + of((©\A)| [E]) + ot (A| [Z]) + ou((0\ A)| [E]) = two. 


Proof. (i) Let am U; ae and ug denote 7* (| [E]), ou* (A| [E]), 1*(©| [E]) and pu*(O| [E]), respectively. 
For every a@ € a 1) the a-cut of pl is 
(3), Gi 


Yr rDE GOT Wy DEy we 


~ 


(14) 


“al(H;| [E]) 


Swapping i with wg we obtain the expression for the a-cut of gu. 
Let X bea fuzzy subset, and let X(‘ and X? be the left and right endpoints of its a-cut. Then 


Keitey =[(%), /(%),- (4), /@)qI 0s) 
because the summations are over empty ranges. On the other hand, for every fuzzy subset X and Va € (0, 1], 
X§ < Xf, and then Va € (0, 1], 1 € “1(O| [E]); ie. 1 € coregl(O| [E]). We can prove that 1 € corequ(O| [£]) 
in a similar way. 7 7 

(ii) Let C C O. Let 1& and u% denote G/*(C| [E]) and pu*(C| [E]), respectively. Then 


(Al [E]) + ot((0\A)| [E)) 
_ i* (Al [E) | i*((@\A)| [E]) 

(*(Al [E]) + o@*((0\A)| [E]) —ol*((0\A)| [E]) + 0*(Al [E)) 

= Ey + Wl a + ata) / Eta + Tala + Batya + TATH\ A) 


(i + TylS\4) [E+ Tm\) (16) 


where k is a fuzzy constant. As = < UG we obtain 7 (A [E]) + o1((@\A)| [E]) < one. The another inequality 
in (ii) is proved in a similar way. 


(iii) It is enough to add the following and together 


Py + 2U4UE\4 + ea 


u (Al [E]) + ou((0\A)| [E]) = 


14 + Talay + By ata + TAG 4 


0.5 An Example In Health Care 


Suppose that we have collected sound level data at a place under study. Also we have interviewed population 
about their annoyance response to noise at this place. Assume that some time in the future, we are interested in 
knowing, only from sound level measurements, without the need of new surveys, if people living at this place, 
is annoyed because of acoustic contamination. 
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Table 0.1: Fuzzy t-numbers in describing interval fuzzy and interval probabilities 


TFN5 a b Cc TEN7 a b Cc 
verylow —1/4 0 1/4 zero -1/6 OO 1/6 
low 0 1/4 1/2 very low 0 1/6 1/3 
medium 1/4 1/2 3/4 low 1/6 1/3 1/2 
high 1/2 3/4 1 medium 1/3 1/2 2/8 
very high 3/4 1 5/4 high 1/2 2/3 5/6 


veryhigh 2/3 5/6 1 
one 5/6 1 7/6 


In estimating annoyance level, we use a fuzzy opted quantization into five exhaustive and mutually exclusive 
states: absent, mild, moderate, intense, or severe. Given sound level measurements from a place, we wish to 
classify people’s annoyance response to noise at this place into one of these five situations according to these 
measurements and previous knowledge. 

Let T; and T> denote the triangular fuzzy numbers (a1, b1, ci) and (a2, bz, cz). The basic operations wich we 
use are: 


TT = (ai ta9,b17 60,0150) (17) 
= = (6 =bis= 04) (18) 
LT, = (1/e1, 1/61, 1/a1) (19) 


where ~ denotes approximation [6]. According to Zadeh’s extension principle, the extension of the tfn 


T = (a,b,c) to the fuzzy sets is given by] 


oe _ J 0 ifl<y 
T(A; a, 6, )(y) = { max{A(a+(b—a)y),A(c+(b—c)y)} ify<1 (20) 
Thus, we naturally extend to 
T((E] “it, Pt) = 7 (Ho; 4, b, c), T(E1;4, b, 0)| Nich (21) 


7 TA: a, b, @); T (Eo; a, b, 0)| Xb<x0 
Tv [min{T (Ho; a, b, C), T(E; a, b, ote Nb eee (22) 


We assume that sound level measurements are reported according to an interval with triangular fuzzy 
numbers as endpoints. These “talking” sound level meters use the term-set TFN{’ = {zero, one, ..., ten}. 
Because of the subjective nature of annoyance, evaluation must be carried out using survey techniques such as 
questionnaires. The term-set TFN5 = {very low, low, medium, high, very high} (Cf. Table|o.1) is used by experts 
and by the general population to classify sound level average regarding to a place. Similarly, in describing 
interval fuzzy probabilities for the final conclusions we consider the term-set TFN7 = {zero, very low, low, 
medium, high, very high, one} (Cf. Table [o.1). Sound level average linguistic estimation and probabilities could 
be represented by any ordered term-pair from TFNs and from TFN7, respectively. 

“Data-based” priors for any annoyance level state could be deduced from information contained in samples 
of past experiences in which the same noise was studied. Assume that the prior distribution and the likelihood 
probabilities are given as shown in Tables [0.2] and [o.3| respectively. For example, a zero-valued likelihood 
probability p(very low | severe) is reasonable, because it means that zero is the probability —regarded as 
representing a degree of reasonable belief or confidence rather than a frequency— that an expert classifies as 
very low the sound level average regarding to a place where population’s annoyance is severe (Cf. Table [o.3). 


’ Let U/ and V be two spaces, f : U — Va function, and S € §(U). The image of S by the extension of f is S’ € §(V), constructively defined as: 
Vy EV, S'(y) = [sup{S(x)|c CEUAy = f(x)} if f-l(y) FS; 0 else]. If 9 is T(a; a, b, c) = max(min((a — a)/(b— a), 1, (e— x)/(c— b)), 0) 
then S’ is T(A; a, b,c) given by 20}. 
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Table 0.2: Annoyance priors 


Population’s annoyance p(-) 
absent very low 
mild very low 
moderate high 
intense zero 
severe zero 


Table 0.3: Likelihoods 


Population’s annoyance — very low low medium high very high 
absent very high very low zero zero zero 
mild very low low low very low zero 
moderate zero very low high very low zero 
intense zero zero zero high low 
severe zero zero zero zero one 


We are interesting in the maximum a posteriori (MAP) estimation. For each interval posterior 
[ol (A;| [E)) , ou (H;| [E])] we calculate the mean fuzzy set of its fuzzy endpoints. Baas and Kwakernaak’s 
fuzzy number ranking procedure [/7] allows us to find the MAP. The final step is to find the nearest words 
—belonging to TFV7— to the fuzzy endpoints of the MAP. This is accomplished by using a distance defined 
on triangular fuzzy numbers. Given T| = (a1, 61, c1) and T2 = (ag, b2, c2), we use an unweighted Euclidean 
format [8]: d(T, T2) = \/(a1 — a2)? + (b1 — 62)? + (cr — €2)?. 

For the data showed in Tables|o.1}0.3] and for an interval fuzzy measurement —reported by a “talking” sound 
level meter— of [six, seven], then, with an interval probability of [high, very high] the best alternative is that 


people living at that place is moderate annoyed because of acoustic contamination. 


0.6 Conclusions 


We have centered on how to learn under a Bayesian point of view from imprecise linguistic data (evidence). 
We have assumed that linguistic imprecision lies at the same time in evidence, priors and marginals, and 
accordingly in the posteriors. We represent this linguistic imprecision using interval fuzzy and interval 
probability assignments. Our linguistic imprecision propagation Bayesian mechanism allows us to express 
conclusions in a natural way such as: “In the light of these results we can say that for a sound level of 
approximately six or seven, then, with a probability, approximately high or perhaps very high, the best 
alternative is that people living at that place is moderate annoyed because of acoustic contamination”. 
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